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ABSTRACT 

The formation of stars is a key process in the early universe with far reaching conse- 
quences for further cosmic evolution. While stars forming from truly primordial gas 
are thought to be considerably more massive than our Sun, stars in the universe today 
have typical masses below one solar mass. The physical origin of this transition and the 
conditions under which it occurs are highly debated. There are two competing models, 
one based on metal-line cooling as the primary agent and one based on dust cooling. 
The recent discovery of the extremely metal poor star SDSS J1029151+172927 pro- 
vides a unique opportunity to distinguish between these two models. Based on simple 
thermodynamic considerations we argue that SDSS J1029151+172927 was more likely 
formed as a result of dust continuum cooling rather than cooling by metal lines. 
We conclude that the masses of extremely metal-poor stars are determined by dust- 
induced fragmentation. 
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1 INTRODUCTION 

The earliest generations of stars ended the so-called cosmic 
dark ages and played a key role in the metal enrichment 
and reionization of the Universe, thereby shaping the galax- 
ies we see today (IBromm & Larson 2004 Glover 2005 Ciardi 



|fc Ferrara|2005| |Bromm et al.|20"09 1. It is therefore remark- 



able that there is still very little consensus on the physical 
processes that govern stellar birth in the early universe. In 
particular, there is still a debate over how and when the 
star formation process switches from the high-mass mode 
that is thought to dominate in the very early univer se (|Abel 



et al. |2002l |Bromm et al.|2002| [Yoshida et al.|2006| [O^Shea 



& Norman 2007 1 to the low-mass mode that presides today 
(e.g., |Kroupa|2002||Chabrier|2003||Bastian et al.|2010|>. 

The relatively new field of "stellar archeology" (Beers 



& Christlicb 2005) attempts to address this issue by com- 
piling observational databases of stars with extremely low 
metallicities in the hope that we can start to constrain the 
various theoretical predictions. The recent discovery of the 
extremely meta l poor star SDSS J1029151+172927 by[Caf] 
fau et alT] ( |2011 1 provides a unique opportunity for this task. 
We recap the basic properties of this star in Section [2] In 
Section [3] we review the two main competing models pro- 
posed to lie behind the transition from the primordial to the 
present-day initial mass function (IMF). These are fragmen- 
tation due to metal-line cooling and due to dust continuum 



cooling. We present simple thermodynamic arguments in 
Section[4]and argue that SDSS J1029151+172927 was more 
likely formed as a result of dust continuum cooling, rather 
than cooling by metal lines. Finally, we conclude in Section 

m 



2 PROPERTIES OF SDSS J1029151+172927 

SDSS J1029151+172927 is a Galactic halo star at the posi- 
tion a = 10 ,l 29 m 15.15 s and S = +17°29'28" (epoch 2000) 
with a g-band magnitude of 16.92 m and an estimated mass 
of about 0.65 Mq (Caffau et al., in preparation). It was se- 
lected from the Sloan Digital Sky Survey and spectral infor- 
mation was obtained with the X-Shooter and UVES instru- 
ments at the Very Large Telescope of the European Southern 
Observatory in Chile Pj 

The star is remarkable in the sense that it is the first 
star to be discovered with elemental abundances in the range 
< 10~ 5 to 10 -4 of the solar value for all of the elements 
measured in its spectrum. This sets it apart from other ex- 
tremely metal-poor stars with [Fe/H]~ —5 which typically 



1 SDSS: www.sdss.org/, X-Shooter: www.eso.org/sci/- 
f acilities/paranal/instruments/xshooter/, UVES: 
www . eso . org/sci/f acilities/paranal/instruments/uves/. 
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have enhanced CNO abundances, such that [C/H], [N/H], 
and [O/H] are all at or above values of roughly —3.5 (e.g. 
Christlieb et ,al.|2004||Besseil et al.|2004[|Nor"ris et al.|2007| 
Frebel et aLpOO Slp] 



3 HIGH-REDSHIFT STAR FORMATION 
3.1 Metal- free Stars 

The first generation of stars, the so-called population III (or 
Pop III) build up from truly metal-free primordial gas. They 
have long been thought to live short, solitary lives, with 
only one extremely massive star with typically ~ 100 Mq 



or more forming in each dark matter minihalo (Abel et al. 
2002}|Bromm et al.|2002l|Yoshida et al.|2006||O'Shea fc Nor- 



man 2007). However, in more recent simulations Turk et al 



(20091 and Stacy et al. (20101 have reported fragmentation 



of primordial gas into a wide binary system for some of their 
cases considered. Widespread fragmentation into multiple 
stellar systems and clusters has also been fou nd by |Clark 



et al. 



(2011 



( [20081 |2011b[ ), |Smith et al\ p0lT| and |Greif et al 
Although these calculations suggest an IMF for pri- 
mordial stars that reaches down to masses perhaps as low 
as 0.5 M©, the primary mode of star formation is expected 
to lead to higher-mass objects with typically several tens of 
solar masses, as the emerging mass function appears to be 
flat ( |Clark et al.[2011a| [Greif et al.|20iTj ). This agrees with 
the analysis of abundance patterns of extremely metal-poor 



stars in the Galactic halo (e.g., Beers & Christlieb 20051, 
which are consistent with the enrichment from core col- 
lapse supernovae from stars in the mass range 20 — 40 Mq 
rather than from pair-instability supernovae of progenitors 



of ~ 200 M ( 


Tumlinson et al.|2004||Tumlinson|2007 


Tomi- 


naga et al. 2007 Izutani et al. 2009; Hcgcr & Wooslcy 2010 


Joggerst et al. 


2010 


1. In addition, the metal enrichment of 



the IGM occurs relatively rapidly following the onset of Pop 
III star formation (see e.g. Greif et al.|2010 Maio et al. 2011 



Wise et al.|2012 l, significantly limiting the total number of 



Pop III stars that are able to form. 



3.2 Metal-Enriched Stars 

The second generation of stars, often termed Pop II. 5, form- 
ing from material that has been enriched from the debris 
of the first stars is the topic of this discussion here. There 
are two main models for the physical processes that deter- 
mine the fragmentation behavior of low-metallicity gas and 
consequently set the stellar mass spectrum. One model is 
based on cooling by atomic lines and the other one by dust. 
Both predict that increasing metallicity leads to a transi- 
tion from a primordial mode of star formation, producing 
predominantly higher-mass stars (even though the charac- 
teristic mass of this mode is subject to debate as discussed 



2 Note that |Caffau et aT] ( |201l| report no direct detection of 
oxygen, because there was no suitable line in the spectral range 
considered. They argue, however, that the spectral fit based on 
one- and three-dimensional model atmospheres is most consistent 
when assuming a 'standard' value of [O/Fe] of +0.6, implying that 
[O/H] 4.4. 



above), to the well-known present-day mode of star forma- 
tion, where the IMF peaks at around 0.2 M Q (e.g., [Kroupa| 



2002 Chabrier 2003 Bastian et al. 20101. However, there 



is disagreement on the physical origin of this transition and 
on the critical metallicity /7 cr it where it occurs. The dis- 
covery of the very primitive star SDSS J1029151+172927 
(Caffau et al. 20111 opens up the unique opportunity, to 
distinguish between these two competing models and to de- 
termine which of these is more consistent with the available 
observational data. 

Metal-Line Cooling — The first model focuses on the 
cooling by atomic fine-structure lines from alpha elements 
such as carbon or oxygen at densities around 10 4 cm -3 . This 
leads to Z CI i t ~ W~ 3 Z@, which can be estimated analyti- 
cally by calculating the metallicity required to produce a 
cooling rate equal to the rate of adiabatic compression heat- 
ing for given halo properties ( Bromm & Loeb 2003 . Santoro 
& Shulll 12006b. This lead iFrebel et al.l (|2007b to combine 



the C and O abundance measures and introduce a transi- 
tion discriminant D tr ans = log 10 (l0 [c/H] +0.3 x 10 [o/H1 ). 
They propose that low-mass stars can only form for Dtrans > 
—3.5. This proposition is challenged by the discovery of 
SDSS J1029151+172927, for which D tr a ns = -4.2, if we use 
the upper limit on [C/H] derived by |Caffau et al. (20111 
and also follow them in assuming that [O/Fe] = +0.6 for 
this star. We note that larger Dtrans values are possible 
if the oxygen abundance is larger than assumed. However, 
Dtrans = —3.5 would require [O/Fe] > +1.9, which seems 
highly unlikely given that C, N and other a elements do not 
show signs of significant enhancement. 

Dust-Induced Fragmentation — Another line of reason- 
ing considers dust cooling as the primary agent for frag- 
mentation jOmukai et"aL]|2005| |2010| |Schneider et al.|2002 



2006 



2012) 



In this model, Z cr it is between 10 to 10 Z@ 
where much of the uncertainty in the predicted value comes 
from our poor knowledge of the dust composition and the 
degree of gas-phase depletion (e.g., [Schneider et al.|[2012[ ). 
Detailed numerical simulations based on a tabulated equa- 
tion of state (EOS ) by |Tsuribe fc Omukai| ( |2006| |2008| > and 
|Clark et al.| ( |2008[ ) demonstrated that indeed fragmentation 
occurs at densities around and above n » 10 12 cm -3 and 
at metallicities 10 -5 — 10~ 6 Zq, leading to an IMF with a 
peak below 1 Mq consistent with the present-day IMF. This 



approach has been refined by Dopcke et al. (2011 \ who per- 



formed similar calculations with time-dependent chemistry 
and found equivalent results. 



4 METAL-LINE VS. DUST COOLING: THE 
CASE OF SDSS J1029151 + 172927 

4.1 Basic Theoretical Considerations 

Our current understanding of stellar birth at the present day 
is based on the intricate interplay between the self-gravity of 
the star-forming gas and various opposing physical agents, 
such as pressure gradients, magnetic fields, or turbulence 
(see, e.g., the reviews by Mac Low & Klessen 2004 McKee 



& Ostriker 2007 Ballesteros-Paredes et al. 20071. Despite 



these complexities, the characteristic mass of the stars that 
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form appears to be tightly controlled by the thermodynamic 



state of the gas (e.g., Klessen & Burkert 



2000 



Bonnell|[2006| |Hennebelle fc Chabrier||2009[ ). That means 



Padoan & 



Nordlund|[2002| |Larson||2005[ |Jappsen et al.|J2005| |Clark fc 



that the thermodynamic properties of the star-forming gas, 
i.e. the balance between heating and cooling processes, are 
the key to understanding star formation in the primordial 
universe. 
In 



Li et al. 



(2003) 



a simple numerical experiment 
showed that turbulent gas with an equation of state P oc n 1 
(where P is the pressure and n is the number density of 
the gas) could fragment efficiently and build up a cluster of 
low mass stars provided that 7 < 1. Conversely, star for- 
mation in gas with 7 > 1 is biased towards forming very 
few and very massive objects. |Larson| |l985 , 2005 ) proposed 
that the transition from a cooling-dominated to a heating- 
dominated regime (i.e. the transition from 7 < 1, which 
implies a decreasing temperature with increasing density, 
to 7 > 1, which implies a rising temperature with ris- 
ing density) marks the point at which the gas fragments 
and therefore is responsible for determining a characteris- 
tic mass scale for fragmentation. Based on this proposal, 
Jappsen et al.| ( |2005[ ) performed a suite of numerical sim- 
ulations using piece-wise polytropic equations of state and 
demonstrated that Larson's proposal appears to be correct 
and that the transition from 7 < 1 to 7 > 1 does indeed 
define a characteristic mass scale for fragmentation (see also 
Bonnell et al.||2006|). 



Physically, this behaviour can be understood as a con- 
sequence of fragmentation in turbulent, self-gravitating gas 
occurring primarily in filaments. Filaments with 7 < 1 are 
unstable to fragmentation ( |Larson 1985 2005] ), while those 
with 7 > 1 are stable and do not fragment. |Kawachi fc| 



Hanawa (1998) showed when studying the similarity solu- 



tion for the evolution of self-gravitating gaseous filaments 
that the collapse becomes slower and slower as 7 approaches 
unity from below and eventually comes to a halt at 7 = 1. 
They suggested that the slow-down of the collapse as 7 ap- 
proaches unity will in reality cause a filament to fragment 
into clumps, because the timescale for fragmentation then 
becomes shorter than the timescale for contraction toward 
the axis of the filament. Simulations of star formation in 
present-day molecular clouds indicate that most stars do in- 
deed form within filaments (e.g. |Smith et al.|2011[ ), and this 



idea also has strong observational support (see e.g. Arzou- 



manian et al. 20111. As numerical models of the assembly 



of the first galaxies show that the gas develops a filamen- 



tary morphology (e.g. Wise fc Aberj|2008 Greif et al.||2008 



Bromm et al.|2009| ), it seems reasonable to assume that fil- 
ament fragmentation will dominate in the earliest galaxies 
just as it does in local star-forming regions. 



4.2 Application to Metal-Poor Gas 

The theoretical picture outlined above suggests that in or- 
der to make predictions for the characteristic mass of stars 
formed in low metallicity gas, we need to establish the con- 
ditions in which the effective equation of state changes from 
7 < 1 to 7 > 1. The locations where such 'dips' in the 
temperature-density (T — n) relation occur are illustrated in 
Figure [I] Note that for an ideal gas with a barotropic EOS 
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Figure 1. A schematic diagram showing the relation between 
the expected equilibrium temperature T and the density n for 
different levels of metal enrichment, Z/Zq = 0, 10 — 5 , 10 -3 , and 
1. The 'dip' in the T — n relation at low densities is due to metal- 
line cooling followed by H2 formation heating, while the 'dip' at 
high densities and subsolar metallicities is due to dust cooling 
followed by the transition to the optically thick regime (r > 1, 
solid red line) for continuum radiation. Thin dashed lines indicate 
the lines of constant Jeans mass. The T — n relation shown here 
is representative of the results from detailed studies. The solar 



metallicity line for n < 10 1 



is an approximation proposed 



by Larson (20051, combined with the expected shift to 7 ft: 1.4 at 
a density of 10 1:j cm~ 3 , where the gas become optically thick (e.g. 
|Masunaga et ak]|l998| . Above n ft 10 15 cm - 3 the temperature 
is high enough to dissociate H2 leading to 7 ~ 1.1. The other 
three curves for the lower metallicity and primordial gas reflect 
the one-zone cal culations by |Omukai et~al] l|2005| |2010|l; see also 
|Schneider et al.| ( |2002| |2006| 2012| , or |Dopcke et al.| | |2011| for a 
three-dimensional calculation. 



(where P oc n 7 ) we obtain the relation 7 = dlnT/dlnn + 1 
and can simply read off the effective index 7 from the slope 
of the T — n relation. Two of sets of features are particu- 
larly important for the thermodynamic analysis here. The 
first sequence of 'dips' is caused by metal-line cooling. This 
can be highly effective at low densities, but becomes much 
less effective for densities n > n cr it, the critical density at 
which the fine structure level populations reach their local 
thermodynamic equilibrium values, leading to an increase 
in the temperature. The growing effectiveness of H2 forma- 
tion heating as we move to higher densities also drives T to 
higher values. As shown by |Omukai et al.| ( [2005] |2010[ ), this 
feature covers the metallicity range < [M/H] < —4 (where 
M denotes the total metal abundance). We see from Figure 
[TJthat minima in the T — n curve induced by metal- line cool- 
ing correspond to Jeans masses significantly above 10 Mq. 



Indeed, Bromm et al. (20011 and Jappsen et al. (20091, who 



directly simulate fragmentation in this regime, report mini- 
mum fragment masses in excess of 10 Mq, with the majority 
of fragments having masses of order 100 Mq . These calcula- 
tions demonstrate that the line-cooling regime is unable to 
directly produce stars with masses around 1 Mq for metal- 
licities of [M/H] ~ —3 and below (see also Smith et al. 2009 1. 
This is the typical mass range of the extremely metal poor 
stars observed so far (Beers & Christlicb 2005 1. We conclude 
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that metal-line cooling cannot explain the physical charac- 
teristics of these highly interesting objects. 

The second sequence of 'dips' in Figure [T] is caused by 
the coupling between gas and dust and covers a range of 
metallicities from the solar value to about Z w 10~ 5 Zq or 
possibly lower if different dust models are considered. It oc- 
curs at densities just below the critical value for the gas to 
become optically thick for continuum radiation as indicated 
by the red line. It is important to note that because this 
sequence of 'dips' lies at much higher densities, the corre- 
sponding characteristic mass for fragmentation falls signifi- 
cantly below 1 M@. Based on this simple thermodynamic ar- 
gument, |Omukai et al.| (|2005||2010|),|Schneider et al.| ([2006 



2012[), and others concluded that the observed extremely 



metal-poor stars in the Galactic halo must be caused by 
dust-induced fragmentation rather than by metal-line cool- 
ing. This proposition was confirmed in numerical simulations 
by|Tsuribe fc Omukai| ( |2006l[2008l ), |Clark et aT] p008) , and 



Dopcke et al. (20111. These calculations show that dust- 



cooling at high densities frequently leads to objects in the 
sub-solar mass range down to metallicities of [M/H] ~ —5. 
This is exactly the physical regime represented by the re- 
cently discovered halo star SDSS J1029151 + 172927. 

A potential objection against this model is that it is 
possible that dust production in very metal poor gas is very 
inefficient (see e.g. Silvia, Smith & Shull 2010), although 



this is still quite uncertain. Nevertheless, if this were the 
case, then the high-density 'dip' would disappear, and the 
dynamical behavior of the gas would be very similar to the 
zero-metallicity case, because metals in the gas phase con- 
tribute little to the overall cooling r ate for [M/H] < —4 
and densities above n 



10 lu cm" J (Omukai et al. 



2005 



|2010[ ). As discussed in Section pj.l| even zero-metallicity gas 
is able to form low-mass stars via disk fragmentation at 



high densities (Clark et al. 2011b Greif et al. 20111. Con- 



sequently, SDSS J1029151+172927 could in principle have 
formed without dust cooling. However, this formation chan- 
nel is thought to be rare and the additional cooling provided 
by dust is needed to build up sufficiently large numbers of 
low-mass stars for the mass spectrum to be consistent with 
the present-day IMF with its steep negative slope from small 
to large masses (Dopcke et al., in preparation). Therefore it 
would provide additional support for the importance of dust 
fragmentation for star formation in the early universe if stars 
like SDSS J1029151+172927 were found to be common in 
our Milky Way. The current attempts to survey larger num- 
bers of very metal poor stars will be able to address this 
issue. 



5 CONCLUSIONS 

In this Letter we have argued that the physical characteris- 
tics (mass and metallicity) of the extremely metal-poor halo 
star SDSS J1029151+172927 are consistent with it having 
formed via dust-induced fragmentation in the high redshift 
universe, but are inconsistent with the proposal of fragmen- 
tation induced by metal-line cooling. 



The fact that Caffau et al. (20111 find one extremely 



metal poor star with [M/H] > —5 out of a sample of only 
six candidate objects suggests that this type of star may not 



be as uncommon as previously thought. With about 2,900 
additional candidate stars in the 7 th SDSS data release, the 
chances of finding more such metal-poor stars seems to be 
quite high. 

A number of previous studies have attempted to com- 
pute the number of low metallicity, low mass stars that one 
would expect to find in the Galactic halo at the present day 



(see e.g. Scannapieco et al. |2006 Salvadori, Schneider & 



Ferrara||2007 Salvadori et al. 12010 1. To derive a very crude 



estimate, we note that of order 10 minihalos with metallici- 
ties below 10 Zq contributed to the early assembly history 



of the Milky Way (Diemand et al. 2005 Gao et al. 20101, 



and that the calculations by Clark et al. ( 2008 1 and Dopcke 



et al. (2011 1 typically produce ten or more low-mass proto- 
stars with masses below 0.8 Mq (the limit for stars to survive 



until the present day; see e.g. Picardi et al.|2 004). Although 
these calculations have a number of serious limitations, if we 
adopt a value of ten low mass stars per halo as a first ap- 
proximation, then we estimate that there are of order of 10° 
primitive stars in our Galaxy. It is interesting to note that 
the value that we obtain from this simple estimate agrees 
to within a factor of ten with the values predicted by more 
complex models (e.g. Salvadori, Schneider fc Ferrara|2007 1. 

Finally, we note that current models of hierarchical 
galaxy formation predict that the first and most metal-poor 
halos to merge will become part of the bulge component of 
the resulting spiral galaxy (e.g., Tumlinson|2010 1, implying 
that this may be the most promising location to search for 
extremely metal-poor and metal-free stars (e.g., White & 
Springcl 2000). In this context it is interesting to note that 
the high eccentricity of SDSS J1029151 + 172927 is consistent 
with it being ejected from the bulge of the Milky Way (Caf- 
fau et al., in preparation). 
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